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Abstract 

If the Higgs-like particle at 126 GeV is the standard model Higgs boson, then 
SM4, the simplest four generation extension of the SM, is inconsistent with the most 
recent LHC data. However, 4G variations (BSM4) are possible if the new particle is 
not the SM Higgs boson and/or if other new quanta modify its production and decay 
rates. Since LHC searches have pushed 4G quarks to high mass and strong coupling 
where perturbation theory eventually fails, we examine the leading nondecoupling EW 
(electroweak) corrections at two loop order to estimate the domain of validity for 
perturbation theory. We find that the two loop hypercharge correction, which has 
not been included in previous EW fits of 4G models, makes the largest quark sector 
contribution to the rho parameter, much larger even than the nominally leading one 
loop term. Because it is large and negative, it has a big effect on the EW fits. It does 
not invalidate perturbation theory since it only first appears at two loop order and is 
large because it does not vanish for equal quark doublet masses, unlike the one loop 
term. We estimate that perturbation theory is useful for mg ~ 600 GeV but begins to 
become marginal for tuq ^ 900 GeV. The results apply directly to BSM4 models that 
retain the SM Higgs sector but must be re-evaluated for non-SM Higgs sectors. 



1. Introduction 

If the Higgs-like particle at 126 GeV is actually the standard model Higgs boson, then 
SM4, the simplest four generation extension of the standard model, is disfavored, first just by 
its small mass which exacerbates the stability [1] and little hierarchy fine-tuning [2] problems 
and second by the most recent production and decay data[3]. However, the virtues of a 
fourth generation remain [1] and the problems of SM4 might be addressed in BSM4 models 
that introduce additional new quanta and/or non-SM Higgs sectors. For instance, although 
rrih *5 vtlq is necessary for vaccuum stability in SM4, two Higgs doublet 4G models can have 
stable vacua with a light Higgs, rrih < rriQ, even for masses as light as ~ 100 GeVpQ. 
Two Higgs doublet 4G models can also be consistent with current LHC data for the 126 
GeV state[Hl El E] and with EW data[E]. In addition, other new quanta could change loop 
mediated Higgs production and decay amplitudes and ameliorate the stability and little 
hierarchy problems. 

The fourth generation can also play a role in scenarios in which the Higgs-like particle 
is not elementary. If it is a strongly bound composite state[9], fourth generation TeV scale 
fermions could be the substrate matter fields on which the new strong dynamics acts. For 
instance, models have been constructed in which conformally invariant strong dynamics act- 
ing on the fourth generation breaks electroweak symmetry |10j. engendering two composite 
Higgs doublets that mix with an elementary doublet [TT] . Similarly, if it is a pseudo-dilaton 
generated by approximately conformal high energy dynamics[12j, a TeV scale fourth gener- 
ation could again provide the matter field substrate. In this case the true Higgs boson must 
be heavy, rrih > 600 GeV, and the 4G fermions provide the necessary oblique corrections to 
ensure consistency with the EW data [13], [lj [Hj. The psuedo-dilaton and Higgs boson could 
mix, as in a similar 4G scenario with Higgs-radion mixing [15]. 

There is then no no-go theorem that definitively excludes a sequential fourth genera- 
tion once our horizon expands beyond SM4. The situation is not unlike that of TeV-scale 
supersymmetry, which is also disfavored in its minimal (MSSM) version but is not ruled out 
in a variety of variations. In both cases, direct searches for the associated heavy quanta, 4G 
fermions or SUSY partners, should be pursued to the limits of the LHC's capability. 

LHC searches for the fourth generation T and B quarks have already pushed the 95% 
exclusion limits to m B > 611 GeV[16] and m T > 656 GeV[TT]. assuming B — > tW and 
T — > bW are the dominant decay modes. If the 4G quarks are stable or decay outside the 
detector the limits are even stronger: using the cross section corresponding to the 737 GeV 
lower limit on long-lived stop production [18], we find a 95% lower limit, rriT > 930 GeVQ^] 

1 We obtained this result from Madgraph[19 with a K factor for gg — > QQ of 1 .3 [20 . 

2 Bound state format ion [21] and/or alternative decay scenarios[22j could allow the 4G quarks to evade the 
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These limits exceed the so-called "unitarity limit" at niq ~ 500 GeV for which the leading 
partial wave amplitude for QQ scattering saturates unitarity in tree approximation [23], due 
to the strong Yukawa couplings of the quarks. This raises the question that motivated the 
work presented here: can perturbation theory be used to obtain the precision electroweak 
constraints on the properties of such heavy 4G quarks? In the process of addressing this 
question we became aware of the importance of the two loop hypercharge correction, that is 
relevant and important even within the perturbation theory domain of validity. 

While sometimes referred to as a unitarity limit, the 500 GeV scale is not a limit on the 
allowed masses but is just a landmark for the onset of strong coupling in QQ scattering. It 
does not precisely indicate where perturbation theory fails in other processes, such as in the 
corrections to the EW data. As in QCD near nonperturbative boundaries, reliability can 
only be estimated process by process, from the magnitude of the higher order corrections 
in each case. We do so here by including the leading nondecoupling two loop corrections. 
We find that the expansion is under reasonable control at tuq = 600 GeV, where the two 
loop corrections to p are roughly ~ 10% as large as the one loop terms, growing to ~ 18% 
at nig = 750 GeV and ~ 25% at ttiq = 900 GeV. For leptons tree unitarity is saturated at 
a higher scale, rriL,N — 1 TeV. [23] Since the electroweak fits typically prefer much lighter 
lepton masses, convergence for large lepton masses is a less pressing issue. 

However, even when the perturbation expansion is a useful approximation, for instance, 
at rriQ = 600 GeV, the one and two loop fits yield very different results. This does not signal 
a breakdown of perturbation theory because it is due to the hypercharge correction, which 
breaks the custodial SU(2) even for equal T and B quark masses, and only first occurs at 
two loop order. To assess convergence we need the three loop hypercharge correction, which 
has not been computed. We use a conservative generic estimate of the magnitude of the 
three loop hypercharge correction, which is consistent with the ratio of the one and two loop 
non-hypercharge contributions, to estimate the effect of the uncertainty it generates. 

Although not surprising with hindsight, it at first seems surprising that the two loop 
hypercharge correction is bigger, and even much bigger in the region of parameter space 
preferred by the EW fits, than the nominally leading one loop quark sector contribution. 
Nondecoupling contributions to p, proportional to the square of the heavy fermion masses, 
require breaking of the custodial SU(2) symmetry that preserves p — 1. These arise at one 
loop order from mass splitting within the quark or lepton doublets, as in the well known top 
quark contribution to the p parameter [2H 123] . Since hypercharge breaks custodial SU(2), it 
gives rise at two loop order to a nondecoupling correction to p, proportional to g' 2 rriQ, even 
if the weak doublet fermions are degenerate. As discussed in section 3, the interplay between 
the lepton and quark contributions to the oblique parameters S and T[22] favors suppression 

above limits. 
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of the one loop quark contribution to T and enhancement of the two loop hypercharge 
correction in the parameter region preferred by the fits. Because the hypercharge correction 
to 5p = aT is large and negative, it offsets positive contributions from quark and lepton 
doublet mass splitting and from CKM4 mixingjH]. Because of the correlation with S the 
preferred region has large L — N and small T — B mass splitting. We use the analytic 
expressions obtained by van der Bij and Hoogeveen[2S] for both the 0((m^ — m 2 B )m 2 Q) and 
0(g' 2 mQ) two loop corrections. 

Our fits incorporate a new two loop result[27j for Rt = T(Z — > bb)/T(Z — > hadrons) 
that causes the p- value of the 3G SM fit to fall to 5% (see section 2). For this study we 
assume negligible CKM4 mixing, at or below the few percent level, for which 4G corrections 
are fully captured by S and T. Like the oblique fit, the x 2 minima for the best SM4 fits 
are typically ~ 2 units lower than for SM3. The SM4 and SM3 fits then have comparable 
p- values, since the SM4 fits, being oblique, have effectively two additional degrees of freedom. 
The two loop SM4 fits have lower x 2 minima than the one loop fits, by ~ 1/2 to 3/4 units, 
since the large negative contribution to T from the hypercharge correction allows the two 
loop fit to approach more closely to the x 2 minimum of the oblique fit, as discussed below. 

The more important difference between the one and two loop fits is in the predictions 
for the 4G masses. As an example of how the fits would be used in practice, we imagine a 
scenario in which the masses of the T, B quarks and charged lepton L are known and use 
the EW fit to constrain the mass of the neutrino N. The resulting differences between the 
one and two loop fits are substantial, even for masses for which the perturbation expansion 
is under control. We also show how the constraint on is affected by the uncertainty in 
the magnitude of the two loop corrections as a function of the quark masses. 

In section 2 we summarize the current status of the SM3 and oblique fits, to establish 
baselines for the SM4 fits. In section 3 we present the fits with the leading two loop non- 
decoupling corrections and compare them to the one loop fits. In section 4 we estimate the 
effect on the EW fit of the uncertainty in the perturbation expansion as a function of the 
4G quark masses. Section 5 is a brief discussion of the results. 

2. Standard Model and Oblique Fits 

We use the data set and methods of the Electroweak Working Group [28]. The SM ra- 
diative corrections are computed with ZFITTER[29], including the two loop SM electroweak 
contributions to m\y and sin 2 6^- eff . The largest experimental correlations are included, taken 
from the EWWG. We use the EWWG data set with one exception: we do not include T w , 
the W boson width, because it is much less precise than the other measurements and has no 
discernable impact on the output parameters. The resulting SM and oblique fits with m# 
fixed at 126 GeV are shown in table 1. 
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Experiment 


SM 


Pull 


Oblique 


Pull 


Aa (5) (m z ) 


0.02750 (33) 


0.02739 


0.3 


0.02754 


-0.1 


m t 


173.2 (0.9) 


173.3 


-0.09 


173.3 


-0.1 


a s (m z ) 




0.1186 




0.1180 




S 








0.05 




T 








0.08 




A LR 


0.1513 (21) 


0.1476 


1.8 


0.1473 


1.9 


A 1 

FB 


0.01714 (95) 


0.01633 


0.8 


0.01627 


0.9 


A e ,T 


0.1465 (33) 


0.1476 


-0.3 


0.1473 


-0.3 


A b 

FB 


0.0992 (16) 


0.1034 


-2.7 


0.1033 


-2.5 


A c 
FB 


0.0707 (35) 


0.0739 


-0.9 


0.0738 


-0.9 


Qfb 


0.23240 (120) 


0.23145 


0.8 


0.23149 


0.8 


T z 


2495.2 (23) 


2495.5 


-0.1 


2496.7 


-0.7 


Re 


20.767 (25) 


20.742 


1.0 


20.737 


1.2 


CTh 


41.540 (37) 


41.478 


1.7 


41.481 


1.6 


Rb 


0.21629 (66) 


0.21475 


2.3 


0.21475 


2.3 


R c 


0.1721 (30) 


0.1722 


-0.05 


0.1722 


-0.04 


A b 


0.923 (20) 


0.935 


-0.6 


0.935 


-0.6 


A c 


0.670 (27) 


0.668 


0.08 


0.668 


0.08 


m w 


80.385 (15) 


80.365 


1.3 


80.383 


0.1 


X 2 /dof 




22.2/13 




20.3/11 




CL( X 2 /dof) 




0.05 




0.04 





Table 1: SM and oblique fits for m H = 126 GeV compared to winter 2012 EWWG data[28]. 

The new Rb calculation causes the x 2 °f both fits to increase by ~ 5 units and the 
p-value of the SM fit to fall to 5%. In addition to Rb the major contributor to the \ 2 is, as in 
the past, the conflict between A h FB and Alr-, which was the principal cause of the marginal 
16% p-value of the fit using the previous Rb calculation. This 95% exclusion should be taken 
seriously, since it is not diluted by a "look elsewhere" effect but is a valid statistical indicator 
of the likelihood that the outliers in the fit could have arisen by statistical fluctuations. We 
should then look either to systematic error or new physics as the most likely explanation. 
Systematic error could be theoretical or experimental, and a leading possibility is the use of 
a hadronic Monte Carlo to assess the effect of gluon radiation on the A b FB measurement [30J . 
With slight modifications most new physics scenarios that addressed the A h FB - A LR conflict 
can also incorporate the new Rb result. See [31] for a model addressing the current fit with 
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Figure 1: 95% confidence region for S and T with ran = 125 GeV. The diamond indicates the x 2 
minimum for the oblique fit and the circle denotes the Standard Model. 

references to the earlier literature. 

Table 1 shows that the tensions are not resolved by oblique new physics, since the 
oblique fit has a similar p-value as the SM fit. Figure 1 displays the 95% CL contour in the 
S — T plane, defined with respect to the best oblique fit, which is at S, T — 0.05, 0.08. The 
95% limit is at x 2 = 20.3 + 5.99 = 26.3. These results change very little if the low energy 
data from Moller scattering and atomic parity violation are added to the EWWG data set. 
The best SM4 fits approach the \ 2 an d S,T values of the oblique fit. 

3. SM4 Fits at Two Loop Order 

The SM4 parameter space is five dimensional, with four fermion masses, my, ms, m^, 
rriL, and the mixing angle #34. We marginalize over various combinations of the five SM4 
parameters and over the three usual SM3 parameters, m t , Aa^ 5 \mz), as(rnz), with m# 
fixed at 126 GeV0 If the SM3 parameters were instead fixed at their SM3 best fit values, a 
procedure employed in many fits of BSM models, we would not obtain the true x 2 minimum 
for the SM4 model, which typically occurs at different values of the SM3 parameters than 
the values in the SM3 fit. 

The fits include the nondecoupling, order (Gpmj) 2 two loop corrections to the T pa- 
rameter, which are protected by custodial SU(2) and vanish when the weak doublet partners 
have equal mass. They have been calculated in a variety of different limits; we use the result 

3 Because mz is known to much greater precision than the other SM3 parameters, the x 2 an d p-values 
are the same whether it is marginalized and constrained or just fixed at its experimental central value. 
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of van der Bij and Hoogeveen|26j. We also use their result for the two loop hyper charge 
correction, proportional to g^GpTn 2 , which is important because it does not vanish for equal 
mass partners, since hypercharge breaks the custodial SU(2). It is computed for equal mass 
partners with rriq = (m^ + m 2 B )/2, resulting in a small, nonleading error, of order a sm 2 9w/^ 
times the one loop term. Because the hypercharge correction is large and only begins at two 
loops, the difference between the one loop correction and the total two loop correction is 
not a valid indicator of the convergence of perturbation theory: the two loop results may 
differ substantially from the one loop results even when the perturbation expansion is under 
control. We will estimate the sensitivity of the fits to the uncertainty in the perturbation 
expansion by scaling the two loop corrections by the generically expected uncertainty. 

For the nondecoupling contributions to the S parameter, which only depend logarith- 
mically on the fermion masses, we use the exact one loop expressions from He, Polonsky, 
and Su[32]. 

In figure 2 we compare one and two loop fits. In the fit on the left mx is varied and 
ttin is marginalized while on the right is varied and ttit is marginalized. In both cases 
rriB = 750 and = 200 GeV are fixed with #34 = 0. The neutrino mass is allowed to 
vary to the 46 GeV lower limit that applies for stable or long-lived neutrinos that escape the 
detector]^] As a function of rriT the two loop fit has a lower, broader, and flatter x 2 minimum 
than the one loop fit and both are approximately symmetric in itit — rriB- Neutrino masses 
at the low end of the allowed range are favored, and the lower \ 2 of the two loop distribution 
emerges primarily at small m^- 

Figure 3 displays the 95% CL contour plots in the {rriL—m N ) - (m T - uib) plane. In these 
plots rriT and are marginalized with uil = 200 GeV, 834 = and m# = 600,750,900 
GeV. The one and two loop contours do not overlap, even for m# = 600 GeV, and become 
increasingly separated as the quark mass is increased. The two loop contours are larger than 
the one loop contours, as is apparent in figure 2. 

The SM4 fits reflect an interplay between the quark and lepton contributions to S and 
T. For #34 = the leading SM4 corrections are given entirely by S and T, so that the oblique 
fit shown in table 1 and figure 1 is the limit for how good an SM4 fit can be. The SM4 fits 
then choose 4G masses to yield 5* and T as close as possible to the the best oblique fit, 
S,T = 0.055,0.08. The quark contribution to S, which is rather insensitive to the specific 
values of m T and m B , is large and positive, Stb — 0.16, well above the preferred value. 
It is offset by the lepton contribution to S which includes a negative term proportional to 
— log(mi/m7v), favoring a large mass difference between L and N with tul > tun as seen 
in figure 2. However large (rriL — tun) induces a large positive contribution to T that can 

4 See [53] for a discussion of this scenario. 



6 




Figure 2: x 2 distributions at one (blue dotdash) and two loops (solid red), varying jut while 
marginalizing over mjy (left) and varying while marginalizing over rrix (right), = 750 and 
rriL = 200 GeV are fixed. The horizontal lines indicate the corresponding 90% confidence intervals. 



force S, T outside the 95% contour in figure 1. The fits then favor small T, B mass splitting 
to minimize the quark contribution to T, and the L, N mass difference strikes a balance 
to achieve the most negative possible contribution to S while keeping T from becoming 
too large. The importance of the two loop quark hypercharge correction is then apparent: 
because it makes a large negative contribution to T it allows (m^/mjv) to increase further, 
resulting in fits that more nearly approach the limiting value of the oblique fit. 

These features are visible in table 2, which dissects the quark and lepton contributions 
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Figure 3: 95% CL contour plots computed with two loop (solid red) and one loop (blue dot- 
dash) corrections with uil = 200 GeV, tub = 600,750,900 GeV, marginalizing over mj- and ttin- 
Diamonds indicate the best fits. 
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m B 


600 


750 


900 


750 : m T [95%] 


m T 


608 


756 


859 


841.4 


m N 


79 


59 


56 


205.5 


Sln 


-0.055 


-0.10 


-0.109 


+0.057 


1 L-N 

rp(2) 

L~N 

T oo 

1 LN 


0.25 
0.0019 
-0.00026 


0.334 
0.0018 
-0.00019 


0.35 
0.0017 
-0.00018 


0.00054 
0.0000071 
-0.0015 


^(TOT) 
1 LN 


0.25 


0.336 


0.35 


-0.00094 


Stb 


0.158 


0.16 


0.164 


0.147 


rpW 

1 T-B 


0.0034 


0.0019 


0.090 


0.45 


m(2) 
1 T-B 
T (Y) 
TB 


0.00040 
-0.14 


0.00035 
-0.256 


0.023 
-0.38 


0.099 
-0.30 


^(TOT) 
TB 


-0.14 


-0.254 


-0.27 


+0.25 


Stot 


0.10 


0.06 


0.055 


0.20 


Ttot 


0.11 


0.08 


0.08 


0.25 



Table 2: First three columns: contributions to S and T for the best fits with mj, = 200 and 
tub = 600, 750, 900 GeV, marginalizing over mjv and mx- Fourth column: tul = 200 and tub = 750 
GeV with rriT at its 95% CL upper limit. 

to S and T. For each fit rriL — 200 GeV is fixed. The first three columns are for the 
best fits with m B = 600,750,900 GeV and the fourth is for m B = 750 GeV with m T at 
its 95% upper limit. For the three best fits the two loop hypercharge correction makes the 
dominant contribution to T, especially for m B = 600 and 750 GeV where the one loop 
quark contribution is negligible in comparison. At the 95% upper limit on my, in the fourth 
column, the hypercharge correction is comparable to although smaller than the one loop 
term. The best fit at m B = 900 GeV reaches the S,T values of the oblique fit. 

4. Convergence of the Perturbation Expansion 

In this section we examine the sensitivity of the fits to uncertainty in the perturbation 
expansion, which increases with increasing quark mass. Because of the large contribution 
from the hypercharge correction, which only begins at two loops, the validity of the expansion 
cannot be judged simply by the difference between the one and two loop fits. 

Consider for instance the fit with m B = 600 GeV in table 2. The one loop result, 
= 0.25, differs by more than a factor two from the total result at two loops , T T ot = 0.11, 
but this difference tells us nothing about the convergence of the loop expansion. The one 
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loop result is completely dominated by the leptonic contribution, for which perturbation 
theory is certainly reliable since ttil and itln are well within the perturbative domain, as is 
explicitly evident from the negligible value of the leptonic two loop contributions to T in the 
table. The quark contribution to T is completely dominated by the quark hypercharge term, 
which is large because it is not suppressed by the near degeneracy of the T and B masses 
favored by the fits. 

To assess the reliability of the expansion we need the next order contribution to the 
quark hypercharge term, a three loop correction that is not known. The best we can do 
is to use a conservative generic estimate of the loop expansion parameter, yq/ATr 2 , where 
Uq = tuq/v is the quark coupling to the Higgs boson and v = 247 GeV. In particular we 
use the known ratio of the one and two loop quark corrections [2"B"] that result from the mass 
difference of T and B for the case \mr — ms\ <C tub, which is consistent with and just a 
factor 3/4 smaller than the generic estimate, 

R _ T T 2 1 B _ 3 tuq 
12 16it 2 v 2 

1 T-B 

We then have R 12 = 0.11, 0.18, 0.25 for m Q = 600, 750, 900 GeV. 

To exhibit the effect of uncertainty of this magnitude we compare fits that shift the 
hypercharge correction by a factor ±i? 12 , 

T m^ r (Y). (1±jRi2) . 

Figure 4 compares x 2 distributions with rescaled as above for ms = 750 GeV and 
rriL = 200 GeV, as in figure 2. The distributions are very similar over most of the allowed 
range in my and mjv, with significant differences only near the upper and lower limits on 
in? and the lower limit on m^. The rescaling of the hypercharge correction shifts the 95% 
limits on my and by 5 GeV or less. 

Similarly, the contour plots in the (m^ — m N ) - (%-mg) plane are compared in figure 
5, for m L = 200 GeV and m B = 600,750,900 GeV, as in figure 3. For m B = 600 GeV 
the three contours overlap quite closely, even though 600 GeV exceeds the tree unitarity 
perturbative limit. For mg = 750 GeV they begin to diverge noticeably, as expected from 
the x 2 distributions in figure 4. Finally, for mg = 900 GeV the three 95% contours are 
almost completely non-overlapping, suggesting that perturbation theory may be of limited 
value at this mass scale. 

Finally we consider an example of how the electroweak corrections will be used in 
practice if evidence of a fourth generation is discovered at the LHC. We imagine that the T 
and B quarks and the charged lepton L are discovered and their masses measured, and we 
consider how well the electroweak fit can then constrain the mass of the yet undiscovered 
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Figure 4: x 2 distributions computed at two loops with rescaled by 1 + R12 (blue dashes) and 
1 — R12 (magenta dotdash), compared with the central value (solid red), mg = 750 and niL = 200 
GeV are fixed. The dotted lines mark the 90% confidence intervals. 



heavy neutrino N. In figure 6 we compare the constraints for one and two loop fits, and 
in figure 7 we compare the two loop fits with Ty B smeared as described above. We assume 
rriL = 200 GeV and mg = 600, 750, 900 GeV. In each case ttlt is fixed at its value for the x 2 
minimum obtained by marginalizing over my and m^, and x 2 is then obtained as a function 
of m N . 

We see in figure 6 that the 90% confidence intervals for of the one and two loop fits 




-200 -100 100 300 -200 -100 100 200 -200 -100 100 200 



m L -m N m L -m N m L -m N 

Figure 5: 95% CL contour plots , marginalizing over my and tun, computed at two loops with 

(Y) 

Tj,g rescaled by 1 + R\2 (blue dashes) and 1 — R\2 (magenta dotdash), compared with the central 
value (solid red). The charged lepton mass is fixed at rag = 200 GeV. Diamonds indicate the best 
fits. 
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m N (GeV) m N (GeV) m N (GeV) 

Figure 6: x 2 distributions as a function of with mi = 200 GeV and my chosen as described 
in the text. The dashed blue and red distributions are from one and two loop fits respectively. The 
intersections of the distributions with the corresponding horizontal lines define the 90% confidence 
intervals for mjv- 



are completely disjoint, even for vtlb = 600 GeV, with increasing separation for larger tub- 
In figure 7 we see that smearing the hypercharge correction by l±i?i 2 results in substantially 
overlapping confidence intervals for itlb = 600 GeV that become almost completely disjoint 
for rriB = 900 GeV. Table 3 shows the effect of the smearing on the 90% confidence intervals. 
For ttib = 600 the impact of the smearing is modest, adding 3 GeV to the upper and lower 
limits. For mg = 900 GeV the effect is appreciable, with the 95% CL upper limit on 
increasing from 65 to 81 GeV. A comparable change would occur for the lower limit but is 
precluded by the direct lower limit onmjy at 46 GeV. 




Figure 7: Two loop x 2 distributions as a function of m^ with mi = 200 GeV and m-p chosen 
as described in the text. The hypercharge correction Tj^ is rescaled by 1 + R\2 (blue dashes) 
and 1 — -R12 (magenta dotdash), compared with the central value (solid red). The horizontal lines 
indicate the 90% confidence intervals. 
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m B 


Unsmeared 


Smeared 


600 


68,90 


65,93 


750 


51,68 


46,76 


900 


48,65 


46,81 



Table 3: 90% confidence intervals for mjy in GeV, as in figure 7. The unsmeared values are from 
the central (red) distributions in figure 7, while the smeared values are the outer envelopes of the 
distributions obtained by rescaling the hypercharge correction as described in the text. 

5. Discussion 

If the Higgs-like particle at 126 GeV is confirmed as the Higgs boson of the standard 
model, then SM4 appears to be excluded, although BSM4 variants with additional new 
quanta might still be viable. If it is a non-SM Higgs boson, e.g., a denizen of a 2HDM, 
then 4G models can also still be viable. We have used SM4 as a laboratory to study the 
effect of two loop corrections on the EW fit for 4G models. The results are qualitatively 
applicable to the broader class of 4G models, but detailed results can only be obtained from 
the examination of each particular model. For instance, for 2HDM models in regions of the 
parameter space with enhanced Yukawa coupling, the effect of the two loop hypercharge 
corrections could be even larger than the already large effect we have found in SM4. 

In the SM4 "laboratory" we have identified an important two loop correction that has a 
big effect on the EW fit in 4G models. Although computed 25 years ago by van der Bij and 
Hoogeveen, it has not previously been included in EW fits of 4G models. It is important 
because it is the first correction that breaks the custodial 577(2) even if weak doublet partners 
have equal mass, and, because it makes a large negative contribution to the rho parameter, 
it allows the fits to approach more closely to the limit of the oblique fit in which S, T are free 
parameters. It is then important in the region of parameter space preferred by the EW fits, 
where the T and B quarks have nearly equal masses so that the one loop quark correction 
is small. In that experimentally preferred region it is by far the largest correction to the rho 
parameter from the 4G quark sector. Because of the hypercharge correction the two loop 
fits differ significantly from the one loop fits even when the quark masses are light enough 
that perturbation theory is reasonably convergent. 

We have also studied the convergence of the perturbation expansion with increasing 
quark mass and its effect on the 4G contraints from the EW data, using a generic estimate 
of the order of magnitude of the three loop hypercharge correction, which is the largest 
unknown term in the parameter region preferred by the fits. We find that for itlq = 600 
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GeV perturbation theory provides useful guidance, even though tree unitarity for elastic QQ 
scattering is saturated at at vtlq = 500 GeV. For ttlq = 900 GeV the estimated uncertainty 
increases so that it undermines the usefulness of the expansion. To obtain better etimates 
of the convergence of perturbation theory, it would be necessary to compute the the three 
loop hypercharge correction in the relevant models. Motivation to meet this challenge could 
be found if (and probably only if) evidence emerges for the existence of a fourth generation. 

This work was supported in part by the Director, Office of Science, Office of High Energy and 
Nuclear Physics, Division of High Energy Physics, of the U.S. Department of Energy under Contract 
DE-AC02-05CH11231 
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